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ABSTRACT 

Using clover improved Wilson fermions we calculate thermal vector meson corre- 
lation functions above the deconfmement phase transition of quenched QCD. At 
temperatures 1.5 T c and 3 T c they are found to differ by less than 15% from that 
of a freely propagating quark anti-quark pair. This puts severe constraints on the 
dilepton production rate and in particular rules out a strong divergence of the dilep- 
ton rate at low energies. The vector spectral function, which has been reconstructed 
using the Maximum Entropy Method, yields an enhancement of the dilepton rate 
over the Born rate of at most a factor two in the energy interval 4 < uj/T < 8 and 
suggests that the spectrum is cut-off at low energies by a thermal mass threshold of 
about (2 - 3) T. 



1 Introduction 



The dilepton spectrum is one of the key observables in the study of thermal proper- 
ties of the hot and dense medium created in heavy ion collisions [0, @, [3|. Thermal 
modifications of the dilepton spectrum have been observed at large invariant masses 
where they reflect the suppression of heavy quark resonances |4| . At low energies it 
is expected that non-perturbative in-medium modifications of the quark anti-quark 
interactions as well as quark anti-quark annihilation in the quark gluon plasma in- 
fluence the thermal dilepton rate. Perturbative 2-loop || calculations as well as 
calculations based on hard thermal loop (HTL) resummed perturbation theory || 
lead to an enhancement of the spectrum at low energies, which dominates over the 
suppression arising from massive quasi-particle poles |], Q or non-perturbative gluon 
condensates [|7| . Indeed, an enhancement of low mass dilepton production has been 
observed in heavy ion experiments ||. 

The thermal dilepton rate describing the production of lepton pairs with energy 
u and total momentum p is related to the Euclidean correlation function of the 
vector current, Jy{r,x) = tp(r, x^^ipij, x), which can be calculated numerically in 
the framework of lattice regularized QCD. This relation is established through the 
vector spectral function, oy(cu,p, T), which directly gives the differential dilepton 
rate in two-flavor QCD, 

dW ha 2 1 
d^ = 27^ ^(e^-l) M ^^ T) ' (1) 

and at the same time is related to the Euclidean vector meson correlation function, 
Gv(t,p,T) = f d 3 x exp(ip ■ x)(Jy(r, x) «/y M (0, 0)), through an integral equation, 

Jo sh(uj/2T) 

where the Euclidean time r is restricted to the interval [0, 1/T]. A direct calcu- 
lation of the differential dilepton rate thus becomes possible, if the above integral 
equation can be inverted to determine ay(uj,p, T). Although finite temperature lat- 
tice calculations, which are performed on lattices with finite temporal extent iV T , 
will generally provide information on GV(r, p, T) only for a discrete and finite set 
of Euclidean times tT = k/N T , k = 1,...,N T , this may be achieved through an 
application of the maximum entropy method (MEM) [[|, |IUf . We will present first 
results on such a calculation in this letter. However, even without invoking sta- 
tistical tools like MEM accurate numerical results on Gy(r,p,T) will themselves 
provide stringent constraints on spectral functions calculated in other perturbative 
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Table 1: Simulation parameters for the calculation of Gv{t,T) on lattices of size 
48 3 x 12 and 64 3 x 16. The third column gives the gauge coupling (3 = Q/g 2 . The 
values used for the clover coefficient c sw , the critical hopping parameter k c and the 
current renormalization constant Zy are discussed in the text. The values for T/T c 
given in the first column are based on results for the non-perturbative /^-function 
given in Ref.(14 . 



or non-perturbative approaches. This in turn will put constraints on the thermal 
dilepton rates. 

We present here results of a calculation of vector current correlation functions 
and reconstructed spectral functions at two temperatures above T c , i.e. T = 1.5 T c 
and 3 T c . We will restrict our discussion to the zero momentum limit (p = 0) of 
these quantities and we will therefore suppress any explicit momentum dependence 
in the arguments of Gy and oy in the following. 



2 Thermal correlation functions 



The correlation functions GV(r, T) have been calculated within the quenched ap- 



proximation of QCD using non-perturbatively improved clover fermions [12, 13 
This eliminates the 0(a) discretization errors in the fermion sector. Moreover, in 
order to explicitly control the cut-off dependence of the numerical results we have 
used lattices of size x N T with fixed ratio of spatial to temporal lattice extent, 
N a /N T = 4, and two different temporal lattice sizes, N T = 12 and 16, respectively. 
Due to the large temporal extent of the lattice, the lattice spacing is rather small 
already at the deconfinement transition and becomes even smaller for larger tem- 



peratures. Using T c = 265 MeV [fTl] to set a scale for the lattice spacing one finds, 



a = 1/N T T ~ (0.75T C /T N T ) fm. For the two different temperatures and lattice sizes 
used here it ranges from a = 0.015 fm to a = 0.04 fm. The comparison of results 
obtained on lattices with different temporal extent and fixed T indeed confirms that 
the results for Gy(r,T) are not significantly influenced by finite cut-off effects. 

The thermal correlation functions have been constructed from quark propagators 
obtained from the inversion of the clover improved Wilson fermion matrix. As we 



3 




1.16 
1.14 
1.12 
1.10 
1.08 
1.06 
1.04 
1.02 
1.00 



G v (xT)/Gv ee ' la '(i:T) 


64^x16 1.5T C 
64 j xl6 3.0T C 


■ 














• 


(b) 




xT 



0.2 0.4 0.6 0. 



0.1 0.2 0.3 0.4 0.5 



Figure 1: Thermal vector current correlation functions, Gy, versus Euclidean time 
r (a) and the ratio of Gy and the free correlator Qfj ee ' lat (b). Shown are results for 
T/T c = 1.5 and 3. The solid line in (a) shows the correlation function for a freely 
propagating quark anti-quark pair, G v ree,lat , calculated with Wilson fermions on a 
finite lattice with spatial extent N a = 64 and temporal extent N T = 16. 



have performed calculations in the high temperature, chirally symmetric phase of 
QCD there are no massless Goldstone modes which could interfere with the con- 
vergence of the iterative solvers for this matrix inversion. We thus could perform 
calculations directly in the limit of vanishing quark masses, i.e. directly at the criti- 
cal values of the hopping parameter k c given in Tab. [I]. The values for c sw as well as 



critical k values were obtained from [|Ti| and interpolation. Note that a finite tem- 
perature critical k defined by a vanishing quark mass might differ from the T = 
value k c by finite lattice spacing and finite volume effects. Indeed, the values for the 



quark mass obtained from the PCAC relation ||T5[ are very small but not exactly 
zero, see Tab. [I| However, since the correlation functions above T c are almost quark 
mass independent close to the chiral limit |HJ it is not crucial to hit k c (T ^ 0) 
precisely. 

Some care has to be taken over the proper renormalization of the vector current 
Jvfi used in the lattice calculations. In order to be able to compare results obtained 
with different iV T and T as well as with calculations performed in different continuum 
regularization schemes we use the renormalized local current j^ orm = (2 k Zy)Jy ll . 
Here Zy(g 2 ) is the current renormalization constant for which we use the non- 
perturbative values determined in |T\|. The overall error on the normalization [T3|, 



O , including the effect of the quark mass not being exactly zero has been estimated 
to be less than 1 %. 

In the high temperature limit the ratio Gy/T 3 can be calculated perturbatively. 
The leading order result is obtained by using free, massless quark propagators in 
the calculation of the current- current correlation function. For vanishing momentum 
of the quark anti-quark pair (p = 0) the corresponding spectral function is given 
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by ay(u,T) = 0.75 tt 2 uj 2 tanh(u;/4T). In this limit the differential dilepton rate 
reduces to the Born rate, 



dW Boin _ 5a 2 1 

dcutfp ^ = 0) = 36vr 4 (e-/ 2T + l) 2 ' (3) 

In Fig. |l]a we show our numerical results for Gy/T 3 calculated at T/T c = 1.5 and 
3. These results are compared to the lattice version of the free vector correlation 
function, G v ee,lat /T 3 ||18|| . It is evident that the correlators obtained at two different 



temperatures agree with each other quite well and stay close to the leading order per- 
turbative (free) vector correlation function already quite close to T c . The deviations 
from Cry 66 ' at /T 3 are better visible in the ratio G y/G v ° e,Ut which is shown in Fig. |l]b. 
It is noteworthy that Gy(r)/T 3 is definitely larger than the corresponding free cor- 
relator, (7 v ee > Iat /j" 3 f or both temperatures and all Euclidean times < rT < 1. This 
is furthermore supported by a calculation of these ratios for different values of the 
cut-off, e.g. for N T = 12 and 16 at the same value of T/T c . Results agree within 
statistical errors. 

As is obvious from Eq. ^| the correlation function Gy(r,T) receives contribu- 
tions from all energies u. Due to the presence of the integration kernel K(t,uj) = 
ch(u;(r — l/2T))/sh(u;/2T) large energies are, however, exponentially suppressed. 
This suppression becomes more efficient with increasing Euclidean time tT. For 
tT ~ 1/2 the correlation function thus is most sensitive to the low energy part of 
ay(uj,T). In fact, using the free spectral function in Eq. ^ we estimate that ener- 
gies larger than to/T ~ 15 contribute only 2% to the central value Gy(l/2T,T)/T 3 
whereas in the low energy regime, uo/T < 3, the corresponding contribution still 
amounts to 12%. Gy(l/2T,T)/T 3 thus is sensitive to modifications of the vector 
spectral functions up to uo/T ~ 15 and provides a stringent test of models and ap- 
proximate calculations of the low energy part of the vector spectral function. We 
find 

GV(1/2T, T) _ j 2.23 ± 0.05 , T/T c = 1.5 , . 

5 V*J 



T 3 \ 2.21 ±0.05 ,T/T C = 3 

which is about 10% larger than the free, infinite volume value Gy ee (l/2T, T) = 2 as 
well as the corresponding free value on a finite 64 3 x 16 lattice, G^ 66 ' (1/2T, T) = 
2.0368. Our result thus suggests that at least for a significant range of energies the 
vector spectral function is enhanced over the free case. Such an increase cannot easily 
be incorporated in simple quasi-particle pictures, which only include massive quasi- 
particle poles in the quark spectral functions. They lead to a reduction of Gy{r,T) 
relative to G v cc (r, T) and in turn to a reduced differential dilepton rate |], 0, |i~9f . 
On the other hand, HTL-resummed perturbative calculations, which incorporate 
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not only poles at non-zero quasi-particle masses but also additional contributions 
from cuts in the HTL-resummed quark spectral functions, show that both non- 
perturbative features of quarks propagating in a hot plasma play an important role 
for the low energy part of the meson spectral functions 0. The cut contributions in 
the HTL approach || as well as in 2-loop perturbative calculations || even lead to 
divergent vector spectral functions at small energies which also makes Gy(r,T) to 
diverge for all Euclidean times r. Assuming that at our aspect ratio of TV 1 / 3 = 4 
finite volume effects do not play an important role, the mere fact that we do find a 
finite result for the vector correlator via Eq. ^implies that ay(u,T) vanishes in the 
limit to — > 0. 



3 Reconstructed spectral functions 

The correlation functions Gy shown in Fig. [1] clearly stay close to the leading order 
perturbative result. This suggests that also the spectral functions ay are close 
to that of the free case. We have reconstructed ay from Gy using MEM, which 



has successfully been applied at T = [pi ITT] and also has been tested at finite 



temperature PH| |. In order to take into account finite lattice cut-off effects in the 



T ) 1 



reconstruction of the spectral function we introduce a lattice version, K l (t,uj, iV, 
of the continuum kernel K(r, to) which appears in Eq. |2|. With the lattice kernel the 
spectral function ay for vanishing momentum is defined through 

poo 

G v (r,T)= / dcu a v (uj,T) K l (t,lu,N t ) , (5) 
Jo 

where Kl(t,uj, N t ) is the finite lattice approximation of a kernel appropriate to 
describe the correlation function of a free boson at finite temperature, 

i at \ 2u; N ^-^ exp(—i27mTT) . . 

K l (t,u,N t ) = — > — - — . 6 

T ^ (2N T sxn(nir / N T )) 2 + {uj/Ty 

This lattice kernel differs from the T = version used in || in so far as we have 
explicitly taken into account the finite temporal extent of the lattice. Of course, 
approaches K in the limit of large A^ r . Moreover, for the energy range important 
for the description of the correlation functions close to tT = 1/2, i.e. for uj/T < 15, 
they differ by less than 2% already on lattices with temporal extent iV r = 16. Using 
this lattice kernel we have checked that the MEM analysis of correlation functions 
constructed from given input spectral functions successfully reproduces these input 
functions. In particular, we find that the free vector spectral function for massless 



6 




Figure 2: Reconstructed vector spectral function cy in units of u 2 at zero momen- 
tum (a) and the resulting zero momentum differential dilepton rate (b) at T/T c = 1.5 
(doted line) and 3 (dashed line). The solid lines give the free spectral function (a) 
and the resulting Born rate (b). The insertion in (a) shows the error band on the 
spectral function at 3T C obtained from a jackknife analysis and errors on the average 
value of (Tv(uj, T)/u 2 in four energy bins (see text). 



quarks can be reproduced already by using only 16 equally spaced points in the time 
interval [0,1/T]. 

In Fig. ||]a we show spectral functions reconstructed from the vector current 
correlation functions at T/T c = 1.5 and 3 on the 64 3 x 16 lattices. In the MEM 
analysis we have taken into account energies up to u/T = AN T and we used the 
lowest order zero temperature perturbative result m(u) = 0.75 n~ 2 u 2 as a default 
model. In order to test the statistical significance of our results we have split our 
total data sets in 8 jackknife blocks and performed a MEM analysis on each of these 
blocks. The resulting jackknife error, which is given as a band in the insertion of 
Fig. |2]a, shows that our data sample is large enough to yield statistically significant 
results. We also have checked that a variation of the default model by 20% only leads 
to minor changes in the low energy part of cy, which stay within the error band 
shown in Fig. ^|a. Another question is to what extent the spectral function obtained 
from the MEM analysis is unique. As the MEM analysis is based on a finite data 
set the reconstruction of the spectral function itself is not unique. This uncertainty 
is incorporated in the MEM-error, which is calculated from the covariance matrix 
of the spectral function |2(J for four energy bins in the interval < u/T < 16. The 



resulting error on the average value of cry (a;, T)/u 2 in these bins is also shown in 
Fig. |a. 

The reconstructed spectral functions show a broad enhancement over the free 
spectral function for u/T > 16, i.e. ua > 1. This regime clearly is influenced by 
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lattice cut-off effects. In fact, a similar enhancement is observed in T = spectral 
functions and has been attributed to unphysical contributions of the heavy fermion 
doublers present in the Wilson (and clover) lattice fermion action [ ITJ . As discussed 



above this regime of large energies does not contribute to the central values of the 
vector correlator for tT ~ 1/2. The relevant energy regime is given in the insertion 
of Fig. |2]a. The enhancement observed for GV(1/2T, T) thus is due to the peak 
in <Jy(uj, T)/u 2 found for uj/T ~ (5 — 6). We furthermore note that for both 
temperatures the spectral functions have a sharp cut-off at small energies. They 
drop rapidly below uj/T ~ 5 and vanish below uj/T ~ 3. This is in contrast to 
perturbative calculations of oy(cu,T), which lead to a divergent spectral function 
for small energies. 

In Fig. |||b we show the differential dilepton rate calculated from Eq. [I] and the 
low energy part of ay shown in Fig. ^|a. The comparison with the Born rate shows 
that for all energies uj/T > 4 the difference is less than a factor 2. For energies 
uj/T < 3 the dilepton rate drops rapidly and reflects the sharp cut-off found in 
the reconstructed spectral functions. Of course, with our present analysis, which is 
based on rather small statistics, we cannot rule out the existence of sharp resonances 
like the van Hove singularities present in the HTL-resummed calculations. It also 
may well be that the broad peak found by us will sharpen with increasing statistics 
as it has been observed in related MEM studies of hadron correlation functions at 
zero temperature 0. Nonetheless, any further sharpening of resonance like peaks 
in a certain energy interval of the spectral function has to be compensated by an 
even closer agreement with the free spectral function in other energy intervals due 
to the constraint given by Eq. f|. 



4 Conclusions 



To conclude, we find that already close to the QCD phase transition temperature, 
i.e. for T = 1.5T C and 3T C , the vector correlation function deviates only by less 
than 15% from the leading order perturbative result. This also is reflected in the 
reconstructed spectral function which deviates by less than a factor two from the 
leading order perturbative form for energies uo > 4T. The most pronounced feature 
of the spectral function and the resulting dilepton rate is the presence of a sharp 
cut-off at low energies, to ~ (2 — 3)T. If a threshold of similar magnitude persist also 
closer to T c , there will be no thermal contribution to the dilepton rate at energies 
uo < 2T C during the entire expansion of the hot medium created in a heavy ion 
collision. This is consistent with present findings at SPS energies || and may become 
visible as a threshold in the dilepton rates for long-lived plasma states expected to 
be created at RHIC or LHC energies. 
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